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We study the evolution of linear perturbations in metric f(R) models of gravity and identify a 
potentially observable characteristic scale-dependent pattern in the behavior of cosmological struc- 
tures. While at the background level viable f(R) models must closely mimic ACDM, the differences 
in their prediction for the growth of large scale structures can be sufficiently large to be seen with 
future weak lensing surveys. While working in the Jordan frame, we perform an analytical study of 
the growth of structures in the Einstein frame, demonstrating the equivalence of the dynamics in 
the two frames. We also provide a physical interpretation of the results in terms of the dynamics of 
an effective dark energy fluid with a non-zero shear. We find that the growth of structure in f(R) 
is enhanced, but that there are no small scale instabilities associated with the additional attractive 
"fifth force" . We then briefly consider some recently proposed observational tests of modified gravity 
and their utility for detecting the f(R) pattern of structure growth. 

I. INTRODUCTION 

Modern cosmology, with its many successes, is built on the concept of an expanding universe governed by the rules 
of Einstein's theory of General Relativity (GR). Assuming a homogeneous and isotropic ansatz for the metric, and 
given the material constituents, one can use GR to determine how the universe and its content evolve with time. It 
has gradually become apparent that, within this construction, ordinary particles and dark matter are not sufficient 
to describe what we observe. Current data strongly favor a universe that exponentially inflated at an early epoch 
and, after a period of decelerating expansion dominated by radiation and matter, recently began accelerating again. 
Within GR, a period of accelerated expansion can only be due to a dominating energy component with a negative 
equation of state, the most common example being a static, or slowly rolling, scalar field. Alternatively, one can try 
to avoid postulating new components by modifying GR. 

In 1979 Starobinsky showed that one can achieve a dc Sitter phase in the early universe by replacing the Ricci 
scalar R in the Einstein action with a function R + const • R 2 [1]. Starobinsky's model, which is now recognized as 
the first working model of Inflation, is also important as the first model for which the quantum metric fluctuations 
leading to a nearly scale- invariant initial density spectrum were discovered [2]. More recently, f(R) models have been 
revisited in the context of ongoing cosmic acceleration, starting with [3, 4] and followed by many others. 

The models proposed in [3, 4] achieved late time acceleration by adding inverse powers of R to the action, i.e. R + 
const -i? - " with n > 1. It was later realized [5-8], that these particular types of f(R) models posses instabilities which, 
among other problems, prevent them from having a matter dominated epoch. Eventually, in [9-11] conditions for 
cosmological viability of f(R) theories were formulated and some explicit models satisfying them were constructed [11- 
15]. In particular, it was explicitly demonstrated in [7, 8, 16] that it is possible to "design" f(R) functions with a 
stable matter era to match any desired expansion history of the universe. 

It is well known that f(R) theories can be recast as GR with an additional scalar field conformally coupled to 
all matter [17-20]. Hence, they are a subset of a wider class of scalar-tensor theories with a fixed value of the 
coupling strength 1 . Since this coupling is of order one in Planck units, it would seem that all f(R) models would 
be automatically ruled out for violating the stringent equivalence principle tests (EPT) in the solar system and on 
Earth [21]. However, just like in the case of Chameleon models [22, 23], the f(R) scalar degree of freedom can be 
exponentially suppressed in regions of high matter density. Hence, while the EPT impose strong constraints on f(R), 
they do not automatically rule them out [11, 24]. In Section IIB we show that the EPT constraints only allow f(R) 
models that are practically indistinguishable from the cosmological constant at the background level. This does not 
necessarily mean that there cannot be differences in the dynamics of perturbations. In fact, as already shown in [11], 
and as we also demonstrate in this Paper, f(R) can lead to interesting signatures that can be seen with future weak 
lensing surveys. 



It is perhaps a matter of personal choice whether or not to call f(R) a "modification of gravity", or treat it as a class of models 
somehow separate from other scalar-tensor theories. At least from the model-building point of view, they are as a class of models that 
operationally do not require postulating a dark energy component. 
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Linear perturbations in f(R) models have already been studied in [8, 25-27] using differing analytical and numerical 
techniques. The results were obtained in either Jordan or Einstein frames, and for particular choices of f(R) models. 
In this work, we identify the general features of perturbation dynamics common to all viable f(R) models and connect 
the results obtained in different frames. In addition, we show that one can interpret the f(R) perturbations in terms 
of the dynamics of an effective dark energy fluid with a non-zero shear. We then demonstrate the signature of f(R) 
as would be seen via some recently proposed observational tests aimed at detecting the cosmic shear [28, 29]. 

We start with a review of f(R) theories and the conditions for their viability in Section II. In section III we solve 
the full set of linear perturbation equations in the Jordan frame in the Newtonian gauge and present the numerical 
solutions, as well as an analytical study of the sub-horizon perturbations in both Jordan and Einstein frames, and 
in terms of an effective dark energy fluid with a non-zero shear. We conclude with Section IV, where we discuss the 
prospects of observationally detecting the f(R) pattern of structure growth. 



II. f(R) GRAVITY MODELS 

We focus on gravity theories described by the action 

S = i/ d4x ^~9 [R + f(R)} + J d^x^—gC^g^] , (1) 

where n 2 = 8ttG and f(R) is a general function of the Ricci scalar, R. The matter Lagrangian, £ m , is minimally 
coupled and, therefore, the matter fields, \i-> iau along geodesies of the metric g^ v . The field equations obtained from 
varying the action (1) with respect to g^ v are 

(l + f R ) R„u - \g» v (R + f) + (g^D - V„V„) f R = k 2 T^ , (2) 

where we have defined f R = df/dR. We take the energy-momentum tensor to be that of a perfect fluid, 

T ia , = {p + P)U^U v + Pg^ , (3) 

where is the fluid rest-frame four-velocity, p is the energy density and P is the pressure. P is related to p via 
P = wp, where w is the equation of state parameter (for pressureless matter w — and for radiation w = 1/3). 
When considering the background cosmological evolution, we take the metric to be of the flat Friedmann-Robertson- 
Walkcr form, ds 2 = o 2 (t)(— g?t 2 + dx 2 ), where r is the conformal time and a(r) is the scale factor. The metric g^ v 
is minimally coupled to matter, hence the stress tensor and its conservation law will be the same as in the standard 
GR. In particular, the continuity equation is the usual 

p + m{p + P) = , (4) 

where an overdot denotes differentiation with respect to conformal time and % = a/a. The f(R) term in the gravi- 
tational action leads to extra terms in the Einstein equations, which now become fourth order differential equations 
(in contrast to the second order of the equations in standard GR). In particular, for our cosmological ansatz, the 
Fricdmann equation becomes 

(1 + f R )n 2 + p \f R + Uf R = ^a 2 p (5) 

and the acceleration equation is 

l - (1 + f R )H 2 + a 2! - + \f R = -^a 2 (p + 3P) (6) 

One can interpret the extra terms in the Einstein equations (2) as a contribution of an effective fluid with energy- 
momentum tensor 

n 2 Tf v = f R R^ - x -! 9liv + ( 9flv n - V^V v )f R . (7) 
The effect of this fluid on the cosmological background can be described in terms of an effective density 
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and the equation of state 



H 2 fR - \o?f - \%f R - \f R 



(9) 



This effective fluid approach can facilitate the comparison with Dark Energy (DE) models, as well as the interpretation 
of the effects of modifications to GR. It is, however, important to keep in mind that the extra terms associated with 
modifications to GR have a geometrical origin. 

By definition, f(R) models with the same matter content and the same w e s will have identical background evolution 
histories. In fact, for any given function w e s(a) one can solve the differential equation (9) and find a family of different 
f(R) with the same expansion history but corresponding to different choices of boundary conditions. In particular, 
for t« e ff = —1, eq. (9) leads to 



f R - 2Hf R + 2 



(u - n 2 ) h 



o 



(10) 



One can then set T-L = (n a /3)[p + p c g] and integrate eq. (10) to find all f(R) that have a ACDM background 
expansion history. Analogous equations can be written for any function w c s {a). 

The role of DE in eqs. (5) and (6) is effectively played by the additional terms due to f(R) ^ 0. There exists, 
however, a complementary, and sometimes conceptually simpler, way in which to approach f(R) theories. It is possible 
to conformally transform the metric and bring the gravitational part of the action to the usual Einstein-Hilbcrt form 
of standard GR. The price one pays for this simplification is a non-minimal coupling between matter fields and the 
transformed metric [30-32] , as well as the appearance of a new scalar degree of freedom playing the role of DE and 
evolving under a potential determined by the original form of the f(R) in (1). Assuming the transformation between 
the two frames is always well-defined, the results obtained using either description should be the same. The frame in 
which matter particles fall along the geodesies of the metric, while the form of the gravitational part of the action 
may be modified, is referred to as the Jordan frame. The frame in which the gravitation part of the action is the 
same as in standard GR, while the matter may be non-minimally coupled to gravity, is called the Einstein frame. 
In the present paper we primarily work in the Jordan frame. When discussing the perturbations, we will make the 
connection between the Jordan and Einstein frame descriptions. The description of the Einstein frame version of f(R) 
theories and the mapping between the two frames can be found in Appendix B. 

The emergence of an additional scalar degree of freedom can be seen directly in the Jordan frame. The role of the 
scalar field is played by f R , dubbed scalaron in [1]. Indeed, the trace of eq (2) can be written as 



Rf R )- Y (p-3P) = — 



(11) 



which is a second order equation for a field f R with a canonical kinetic term and an effective potential V c ff(f R ). By 
design, the f(R) theories we consider must have \f <C R\ and \f R \ <C 1 at high curvatures to be consistent with our 
knowledge of the high redshift universe. In this limit, the extremum of the effective potential lies at the GR value 
R = k 2 ( p — 3P) . Whether this extremum is a minimum or a maximum is determined by the second derivative of V e s 
at the extremum, which is also the squared mass of the scalaron: 



d 2 V cS 
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At high curvatures, when \Rf RR \ <C 1 and f R — > 0, 
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It then follows that in order for the scalaron not to be tachyonic one must require f RR > 0. Classically, f RR > is 
required to keep the evolution in the high curvature regime stable against small perturbations [5, 33]. The scalaron 
mediates an attractive "fifth force" , which has a range determined by the Compton wavelength 
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While X c is large at current cosmological densities, terrestrial, solar and galactic tests are not necessarily violated 
because the scalaron acquires a larger mass in regions of high matter density. This is essentially the Chameleon 
mechanism of [22, 23]. 
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A. Designer f(R) 



As mentioned earlier, the 4th order nature of f(R) theories provides enough freedom to reproduce any cosmological 
background history by an appropriate choice of the f(R) function [8, 34, 35]. Instead of using eq.(9), we integrate the 
modified Friedmann equation (5) directly to find f(R) models with specified expansion histories. In this subsection 
we describe this "designer" procedure, which generalizes that of [8] to include a radiation component and to allow for 
a time-dependent w e g(a). 

As in [8], let us introduce dimensionless variables 



f(R) H 2 , N 

where H = a~ 1 da/dt 1 t is the physical time, and H is the Hubble parameter today. We fix the desired expansion 
history to that of a flat universe containing matter, radiation and dark energy with a given equation of state. Namely, 
we set 

E = n rn a- 3 + fl r a~ 4 + p e s/pc = E m + E r + E cS (16) 

with 



E c ff = (1 - fl m - fl r ) exp 



J a 



— 31na + 3 / w e s(a)dln a 



(17) 



H2 - 3(425 + . (18) 



where p° c = 3Hq/k 2 is the critical density today, fl m ,r are present day fractions of matter and radiation densities and 
we have defined Ei = p^/ p° cr (that will be used later in the perturbation equations) . The other background functions 
that appear in the modified Friedmann equation (5) can be expressed in terms of E and its derivatives. For example, 

R_ 

1 a 

where the prime denotes differentiation with respect to In a. Substituting (16) into eq.(5) yields a second order 
equation for y: 

y "-{ 1 + §E + W) V ' + 6^E y = Se- EcS ' (19) 

In order to set the initial conditions, let us consider the general and the particular solutions of (19) at early times, when 
the effects of the effective dark energy on the expansion are negligible. At a certain early value In cii, the homogeneous 
part of (19) is satisfied by a power law ansatz y cx a?. Substituting this ansatz in and solving the quadratic equation 
for p yields 

P± = ^[-b±Vb 2 -4c) , (20) 

(21) 



where 

7 + 8r 



2(1 + n) ' 2(1 + n) ' 

and a eq is the scale factor at the radiation-matter equality. The decaying mode solution corresponding to p- leads to 
a large f(R) at early times which makes it unacceptable, and we set its amplitude to zero. The particular solution at 
In di can be found by substituting 

y p = A p E cS ( ai ) (22) 
into (19), where E e g is given by (17). One then finds 

-6c 

-3«4f + 9w 2 s + (I8~b)w cf[ + 9 - 3b + c ' 

Put together, the initial conditions at lnaj are 

y% = AeP+ laa >+y p (24) 

y[ = P+ AeP+ lna * - 3[1 + w cS ( ai )]y p , 

and A is the remaining arbitrary constant that can be used to parametrize different f(R) models with the same 
expansion history. 
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B. Viable models 

While it is possible to find f(R) models that reproduce any given expansion history, they must satisfy additional 
conditions to be consistent with existing experimental and observational data. In this subsection we summarize these 
requirements and examine their implications for the acceptable range of values of w e s. 

The most common argument against f(R) theories is based on their identification with scalar-tensor gravities [17- 
20]. In particular, it has been shown that in the metric formalism considered here, f(R) corresponds to a Brans-Dickc 
theory with lobd = [20]. Since low values of ujbd have long been ruled out by solar system tests, it would appear 
that any f(R) theory is automatically ruled out. However, it is also known [11, 22, 23] that in regions of high 
concentration of matter, the non-minimally coupled scalar degree of freedom acquires additional mass and, therefore, 
can be strongly suppressed. In such situations, simple arguments based on the formal equivalence with Brans-Dickc 
are insufficient to make conclusive statements about viability of f(R) theories. Instead, one can formulate a set of 
requirements for these theories to satisfy in order to pass all of the tests. While these requirements are quite strict 
and generally discourage large deviations of f(R) from the trivial case of a cosmological constant (/ = — 2A), there 
can be room for non-trivial observable differences. Such differences are interesting to understand in some detail and 
work has already been done on this in [8, 27, 36] 

Below we list and discuss the conditions that f(R) models of cosmic acceleration must satisfy in order to be viable. 

1- fun > for R » fan. Classically, this follows from requiring the existence of a stable high-curvature regime, 
such as the matter dominated universe [5]. Quantum mechanically, mj^ w Jrr > ensures that the scalaron 
is non-tachyonic. 

2. 1 + /r > at all finite R. The most direct interpretation of this condition is that the effective Newton constant, 
G e ff = G/(l + fa), is not allowed to change sign. Among classical implications of 1 + fn < is the universe 
quickly becoming inhomogeneous and anisotropic [37, 38] . The quantum mechanical significance of this condition 
is in preventing the graviton from turning into a ghost [39] . 

3. fa < 0. Given tight constraints from Big Bang nucleosynthesis (BBN) and the Cosmic Microwave Background 
(CMB), we want regular GR to be recovered at early times, i.e. f(R)/R and fa — > as R — > oo. Together with 
fun > 0, this implies that fn must be a negative, monotonically increasing function of R that asymptots to 
from below. 

4. fn must be small at recent epochs. This is not required if the only aim is to build a model of cosmic acceleration, 
without trying to satisfy the solar and galactic scale constraints. This condition ensures a small difference 
between the value of the scalaron field in the high density galactic center and the value in the outskirts of the 
galaxy [11]. As shown in [23], the difference between these values is effectively the potential difference that 
sources the attractive "fifth" force acting on objects in the vicinity of the galaxy. Analysis of [11] suggests that 
the value of \fn\ today should not exceed 10~ 6 . However, as stressed in [11], this bound assumes that in f(R) 
models galaxy formation proceeds similarly to that in GR. Hence, while it is certain that \fu\ must be small, 
any specific bound on its value today will be unreliable until galaxy formation in f(R) is studied in N-body 
numerical simulations. 

These restrictions reduce the space of allowed values of w c g- Namely, there are functions w e s(a) for which none of 
the solutions of eq.(19) satisfy all of the above conditions. To understand this better, let us examine a typical solution 
of eq.(19) for a given expansion history. It can be written as a sum of the solution of the homogeneous part of the 
equation plus the particular solution. During the radiation and matter domination the homogeneous part of eq.(19) 
is only weakly dependent on w e s. In that regime, the main impact of w e s(a) is through the driving term on the RHS 
of eq.(19) which sets the time dependence of the particular solution. To gain intuition into which w e s are viable, one 
could first "turn off" the homogeneous part by setting A = in (25) and see if the above Conditions are satisfied. 
Then, if any of the Conditions are violated, one could try to remedy the situation by a suitable choice of A. 

The particular solution is proportional to p ff, whose absolute value is a decreasing (increasing) function of time 
for models with w e s > — 1 {w e s < —1). In the special case of w e s = — 1, the particular solution is a constant. In this 
case, all positive values of A result in solutions satisfying Conditions 1 and 3, while all solutions with negative A have 
fun < and fx > 0. One can also satisfy Condition 4 by choosing an appropriate A in order to make as small 
as needed. Three sample solutions for w e s = — 1 with A = 0, A < and A > are shown with black solid lines in 
Fig. 1. A = corresponds to the case of the cosmological constant, f(R) = — 2A. 

Let us next consider models with a constant w c g ^ — 1. Particular solutions in the case of w e s < — 1 are unacceptable 
because they lead to f(R) that violate Conditions 1 and 3, as can be seen from the red short-dashed line in Fig. 1. 
One could try to remedy the situation by choosing a sufficiently large positive A, to make Jrr > 0. However, all such 
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FIG. 1: f(R) (top panel) and fn (bottom panel) vs R obtained using the designer approach for several given expansion histories. 
The three solid black lines correspond to w cB = -1 with A = (f=-2A), A < (-//2A > 1, fn > 0) and A > (-//2A < 1, 
/fl < 0). The red short-dashed line corresponds to tu e g = —1.01 with A = 0. The blue dot-dashed line is for w e g = —0.99 with 
A = 0, while the blue dot line corresponds to w e s = —0.99 but with A adjusted to make \/r\ — 10 -6 today. The blue dot on 
the lower panel corresponds to 10 x fn which makes it easier to see. Finally, the green long dash denotes a model with A = 
and a time varying w e s (a) which crosses from —0.99 at early epochs to —1.01 for z < 1. 



attempts lead to models that initially have fun < and then cross through fun — at some finite R. Such models 
would have instabilities and hence w c g = const < — 1 is not allowed. 

Models with w c s — const > — 1 are generally better behaved, since the corresponding particular solution satisfies 
Conditions 1 and 3 as can be seen by looking at the dot-dashed blue line in Fig. 1. One does run into troubles, 
however, when trying to satisfy Condition 4 by making \fn\ small today with a suitable choice of a negative A. Then 
the particular solution will still force \fn\ to grow in time at high R, while the homogeneous part will bring it down 
to a small value today. Such fn will have a turning point and, hence, will violate both Conditions 1 and 3. This case 
is illustrated with a blue dotted line in Fig. 1. However, if one chose to only study cosmological implications of f(R) 
and ignore the galactic and solar constraints, then histories with w e g = const > — 1 would generally be acceptable. 

One gains more flexibility by allowing for a time-dependent equation of state. In particular, it is possible to 
construct f(R) that have w e g crossing —1 without violating Conditions 1 and 3. An example is shown in Fig. 1, 
where the long-dashed green line corresponds to a model that has w e g = —0.99 at early times and transitions to 
-1.01 at redshift of 1. 

Quite generally, Condition 4 discourages large deviations of w D g from —1, since otherwise one must deal with turning 
points in fji at finite values of R. Hence, if we want to have reasonable physics on solar and galactic scales, we must 
have w e g ps — 1. As a consequence, if future data actually finds any deviations from w e g = —1, it would automatically 
rule all f(R) models out. Only if galactic and solar tests are ignored, a richer variety of expansion histories can be 
considered. For the numerical analysis in the rest of this Paper we will only use models with w e f{ = — 1. 

Let us remark that, while such designer models appear extremely contrived, they are not necessarily more fine- 
tuned than quintessence dark energy models in terms of the number of free parameters. One could argue, in addition, 
that a universal order one coupling of the scalar field to matter is better motivated by high energy physics, where 
dilaton fields can appear naturally. Producing minimally coupled scalar fields requires additional fine-tuning. From 
the observational point of view, a quintessence dark energy with w — > — 1 would be impossible to distinguish from A. 
The f(R) models, on the other hand, predict a ACDM background expansion but detectably different dynamics of 
structure formation. 
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An additional potential problem with f(R) models of "dark energy" was recently pointed out by Starobinsky [14]. 
He observed that the perturbation 6R of the Ricci scalar in f(R) models in the high curvature limit is a rapidly 
oscillating function of time, with an amplitude that can easily be larger than the average value of R. An oscillating 
background could lead to an overproduction of scalaron particles that would impact the Big Bang nucleosynthesis. 
Commenting on the quantum particle production would go beyond the scope of this Paper, although this is an 
important issue that should be further examined and clarified. At the classical level, however, it would seem that 
large rapid oscillations of R do not necessarily imply a problem. R is constructed from derivatives of the metric, 
so perturbations of the metric can remain tiny even when their derivatives are large when they undergo very rapid 
oscillations. In fact, as we show in the next Section, the metric potentials and the density perturbations are well 
behaved, despite SR being large and oscillatory. In the equations for the metric, SR always appears pre-multiplied by 
the small factor f RR which keeps the metric perturbations in the linear regime. 



III. PERTURBATION THEORY IN THE JORDAN FRAME 



We consider scalar perturbations in the Newtonian gauge, with the line element given by 

ds 2 = -a 2 (l + 2*)rfr 2 + a 2 (l - 2^>)dx 2 , (25) 

where ^ = ^(t, x) and $ = <I>(t, x) are small time and space-dependent perturbations to the metric. In the rest of 
the paper all the perturbed quantities and the equations are presented in Fourier space, where the different /c-modes 
evolve independently. Primes will indicate derivatives with respect to In a and quantities such as H, k and f RR are 
expressed in units of H (or Hq) so that they are dimensionless. 

For a general fluid component, the energy-momentum tensor has the following expansion at first order 

T ° = -p(l + S) , 

lf = -(p + P)v i = -pV i , 

T i J = {P + 5P)5)+-K), (26) 

where 8 = 5p/p is the density contrast, v is the velocity, V = (1 + w)v, 5P the pressure perturbation, and we define 

the anisotropic stress pH = (k j h — ^Sj^j 7r*-, where ttj denotes the traceless component of the energy-momentum 

tensor. In the Jordan frame, matter is minimally coupled and follows the geodesic of the usual metric g^ v . Thus, 
there is no explicit dependence of the matter Lagrangian £ m [%j, g^] on the f(R) modifications, and the conservation 
equations for matter do not differ from the conservation equations of standard general relativity. We refer the reader 
to Appendix A for the explicit expression of these equations, as well as for the full set of modified linear Einstein 
equations. Here we want to focus on the anisotropy and Poisson equations; the former relates the two Newtonian 
potentials ^ and $, while the latter describes how the curvature potential $ depends on the matter comoving density 
perturbation 

A = 6 + 3^-V. (27) 
k 

For theories described by the action (1), the anisotropy equation reads 

^-^=~E i U i -f R (^-^>) + f RR SR, (28) 

where SR is the linear perturbation (A7) to the Ricci scalar R, repeated indices denote a sum over the matter fields, 
and we have defined Ei in (16). The Poisson equation, obtained combining equations (A3) and (A4), is 



k 2 3 A 

— $ = --E^ 
a z 2 



Ir^ ~ l^f^ SR + \ H2 f'n (* + $ ') + \HH'f RR 5R 



(29) 



Using equation (28), and neglecting any anisotropic contribution from matter fields (i.e. setting = 0), we obtain 
the following relation between the Newtonian potentials 

$ _ ^ = l^-SR (30) 

F 



where we have defined F = 1 + 
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FIG. 2: These plots show the numerical solutions for the transfer function <J> +/ /<E>™ tmi ; the ratio of the potentials and 
the slip x = IrrSR as functions of the scale factor a for different scales for an f(R) model with w e s = — 1 and = — 10 -4 . 
The dotted line corresponds to the scale k = O.Olh/Mpc, the short-dashed line to k — O.lh/Mpc and the long-dashed line to 
k = 0.5h/Mpc. Finally, the solid line shows the scale-independent behavior in the ACDM case. The right panel shows the tiny 
oscillations in \ a t early times, shown for k — O.lh/Mpc. 



In standard GR, once we neglect any matter shear, the anisotropy equation is simply a constraint, >!' = $, which 
reduces the number of independent perturbed variables. The Poisson equation is also an algebraic relation between 
the curvature perturbation $ and the matter density perturbation E^A,;. 

On the contrary, in f(R) gravity, eqs. (29) and (30) are dynamical. Indeed, due to the higher order of the theory, 
the linear equations contain extra dynamics. In the specific case of the anisotropy equation, the extra dynamics is 
encoded in the slip between the gravitational potentials. We will choose this slip as one of the perturbed variables 
to evolve. To render the equations more treatable numerically, and also to facilitate the physical interpretation of 
the results, we introduce new variables, functions of the Newtonian potentials $ and \P. Namely, along with the 
perturbations of matter fields, we choose to evolve the following variables 

$ + \& 

= (31) 
X = f RR SR = F(* - *) . (32) 

Any non-zero value of the variable \ will signal a departure from standard GR. The variable $ + is the combination 
of potentials that affects propagation of light, leading to the Integrated Sachs- Wolfe (ISW) effect in the CMB and 
Weak Lensing (WL) of distant galaxies. 

In terms of the variables (31)-(32), the anisotropy equation (30) becomes a constraint (which allows us to eliminate 
$ in terms of $ + and x). We ignore the radiation component due to its relative unimportance for the late time 
dynamics of perturbations. We can then obtain two coupled first order differential equations for \ an d < £ > + from the 
momentum (A4) and the Poisson (29) equations 



3 aE m V m f l-F'N 3_F 
2 HkF { 2 F J + 4F 2 ' 



$' + = ^^^-(l+~]$ + + 7T5X (33) 



, 2E m A m F ( F H F\ , ( 2 k F \ 

X = + i + 2 x ~ 2-F$'+ -2F\l+ — $+ . (34) 

A H 2 F f \ F H F J + V 3a 2 H 2 F'J + V ' 

The numerical solutions of the above equations are shown in Fig. 2 for an f(R) model with w e g = —1 and f R = 
fn{z — 0) equal to — 10~ 4 . Such a value of f R is outside the bounds of Condition 4 of Subsection II B. However, we 
choose to use this model in all of our plots because the pattern of f{R) modifications can be easily seen "by eye", 
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while all qualitative features of the perturbation dynamics are still the same as in models with \f R \ < 10 -6 , that are 
more likely to satisfy Condition 4. We start evolving the equations at z = 1000, when the deviations from GR are 
small, and set $ + = — 1 and \ — at the initial time. We then use the standard GR relations 

2k 2k 2 
Vm = 3^*+ ' Am = "3^* + (35) 

to set the initial conditions for the matter perturbations A m and v m . 

At early times, when the effect of f(R) on the background expansion is negligible, and when the wavelength 
corresponding to k is larger than Ac (but still inside the horizon), the potential stays constant as one would expect in 
GR. As soon as the mode enters inside the Compton radius, the additional attractive force mediated by the scalaron 
starts to enhance the growth. The potential eventually starts to decay when the universe begins to accelerate. The 
plot of $+ in the left panel of Fig. 2 clearly shows these two trends - the enhancement due to the "fifth force" and 
decay due to an accelerating background. In the right panel of Fig. 2 we show a magnified plot of the evolution of 
the same quantities at early times. One can see that \ oscillates rapidly, but the amplitude is very small. These 
oscillations correspond to very large oscillations of SR = x/fun, (also pointed out in [14]), but they do not have a 
noticeable effect on the evolution of the metric perturbations. 

To better understand the asymptotic values of \ and in Fig. 2, we analytically examine the evolution of 

perturbations in the Jordan and Einstein frames in the sub-horizon regime. We then interpret the modified dynamics 
of perturbations in f(R) in terms of an effective dark energy fluid with a non-zero shear. 



A. The behavior of perturbations on sub-horizon scales 



We arc interested in the behavior of matter and metric perturbations for modes that are inside the horizon well 
after the radiation-matter equality. Therefore, in the following analytical study we neglect radiation and deal just 
with cold dark matter (CDM) characterized by w = 0. It is a well known result of standard GR that the CDM 
density perturbations grow linearly with the scale factor on sub-horizon scales during matter domination, while the 
potentials $ and * stay constant. When the cosmological constant begins to dominate, the growth of structures slows 
down and the Newtonian potentials start decaying, resulting in a late-time Integrated Sachs- Wolfe (ISW) signal which 
contributes to the low multipoles of the CMB. These results can be easily derived by considering the sub-horizon, 
fc/a» H, version of the evolution equations. 

In f(R) theories, however, there is an extra length scale associated with the additional scalar dynamical degree 
of freedom - the Compton wavelength Ac (14) of the scalaron. This scale separates two regimes of sub-horizon 
gravitational dynamics [11, 40], so the sub-horizon approximation must be taken with more care. On scales A » Ac, 
the scalar field is massive, the "fifth force" is exponentially suppressed, thus deviations from GR are negligible. 
However, on scales below the Compton wavelength, the scalaron is light and deviations are significant. The relations 
between $ and \I>, and between them and the matter density contrast will be different below the Compton scale and 
that affects the growth rate of structures. Different approaches to the study of the growth of structures have already 
been employed in the literature [8, 41, 42]. It is instructive to compare these approaches. We do it in the following 
subsections, where we also introduce a new method in the variables \ an d , as well as perform a parallel analysis 
in the Einstein frame. 



1. In Jordan frame 

On sub-horizon scales, the equations for $ + and \ reduce to 

k 2 , 3 E m 8 m 



- (36) 



2f RR k 2 
3f RR k 2 + a 2 F' 



X - ~ J ™, _*„ F* + , (37) 



where we have used A ~ 6 on sub- horizon scales. From (37), it is clear that the behavior of \ is characterized by a 
transition scale set by the Compton wavelength of the scalaron Ac (14). 

It is useful to introduce a parameter Q, approximately defined as the squared ratio of the Compton wavelength to 
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FIG. 3: The evolution of the growth factor for the CDM [A m (fe, a)/a]/[A m (k, «i)/ai] as a function of redshift z and scale k. 
The left panel corresponds to an f(R) model with w e s = — 1 and f% — — 10~ 4 . In the right panel we show the corresponding 
ACDM pattern for comparison. One can see the scale-dependent behavior of the growth factor in f(R) as opposed to the scale- 
independence of the ACDM case. The dashed line crossing diagonally on the left plot corresponds to the Compton transition 
wavelength given by Q — 1 (38). Deviations from ACDM become important on scales below Ac- For modes with A < Ac 
during matter domination, there is an enhancement in the growth due to the "fifth force" introduced by the modifications to 
GR. Eventually, the universe starts accelerating and the growth slows down. However, in comparison to the ACDM case, such 
slowing is delayed in a scale dependent way. 



the physical wavelength of a mode 2 



Then we can rewrite eq. (37) as 



3*!/rk r /Ao\ (38) 



F V A 



*HtV* + - (39) 

For scales A ^ Ac, the parameter Q — > and the variable x is negligibly small. However, for scales below Ac, Q 3> 1 
and x becomes of the same order as <I> and asymptotically approaching x — —2/3F&+. This can also be seen from 
the plots in Fig. 2. 

Next, we want to understand the impact of the relation (39) on the evolution of the gravitational potentials and 
the growth of density perturbations. We choose to do it in the variables <£> and to render the comparison with other 
literature easier. From eq.(39) it follows that 

* * ^* . (40) 
3 + 2Q v ' 

On scales A 3> Ac, Q — > 0, and the standard relation \P ~ $ still holds. However, on scales below Ac, the relation 
between the metric potentials becomes \& ~ 2$. 

In order to see how the growth of structures is affected by f(R) modifications, let us combine eqs.(Al)-(A2) into a 
second order differential equation for the matter density contrast. On sub-horizon scales, considering only CDM, the 
equation is 

C+(l + f^K + ^* = 0. (41) 



2 A similar parameter Q was first introduced in [41]. Here we choose to have it positive and to correspond to the ratio of the wavelengths. 
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FIG. 4: The function probed by the cross-correlation of large scale structure with ISW, A m ■ d&+/dz, as a function of scale k 
and redshift z. The left panel corresponds to an f(R) model with w c s = —1 and /jj = — 1CP 4 and shows a characteristic scale- 
dependent pattern. The right panel corresponds to ACDM. The dashed line crossing through the left panel corresponds to Q — 1 
(38), i.e. it corresponds to the Compton wavelength of the scalaron Ac (14). For scales A < Ac, during matter domination, one 
can clearly notice the effect of the "fifth force" which suppresses the cross correlation and can actually make the correlation 
negative. Therefore, a negative cross correlation signal at early redshifts (corresponding to matter era), is a signature of f(R)- 
The acceleration of the background will eventually contrast the "fifth force" and lead to a positive cross-correlation. 



Combining the sub-horizon version of the Poisson equation (29), with (40), we can derive an expression for the 
potential VP in terms of the matter density perturbation 

3 1_13 + 4Q 

a 2 2F31 + Q V ' 

One can interpret the modification as the Newton's constant being rescaled in a time, scale and model dependent way, 
inducing a scale (and model) dependence in the growth of structures. In particular, for a mode inside the horizon 
but well above the Compton wavelength, Q <C 1, the only effect of the f{R) modifications to the growth of structure 
is a 1/F reseating of the constant G. But when the mode is well inside the Compton scale, the Newton's constant is 
rescaled by a factor 4/(3F). Hence, we have two regimes separated by the scalaron Compton wavelength: 

G 

~ $ , X - , G off ~ — , for A > A c 

*~2$, X ~-~F$ + , GeffSi— , for A « A c (43) 

This effect can be clearly seen in Fig. 3. On scales smaller than the scalaron Compton wavelength, the modifications 
introduce a fifth force which would enhance the growth of structures. The rate of growth will depend on the balance 
between the fifth force and the background expansion. For modes that cross below Ac during matter domination, 
the effect of modifications is maximized, as the potential $ + grows in the absence of background acceleration. When, 
however, the background expansion starts accelerating, it compensates for the enhancement due to the fifth force. 
Eventually, potentials start to decay, but at a lesser rate than in the ACDM model. 

A characteristic observational signature of f(R) models would be an ISW effect during the matter era if one could 
correlate the distribution of large scale structure at z > 2 with the CMB. In fact, such a correlation would be negative, 
since potentials would be growing, and not being constant as in the usual case. This is clearly seen in Fig. 4. One 
must be realistic, however, and keep in mind that large statistical errors associated with the ISW measurements, and 
the smallness of the effect expected in viable f(R) models, will likely render the ISW based tests of f(R) useless. 
Unless, of course, missions such as LSST [43] and SKA [44], which will be able to probe structures at redshifts z < 3 
and z > 3 respectively, find a statistically significant positive LSS-ISW correlation signal at high redshifts, which 
would effectively rule out f(R) (and ACDM). 

From eq.(36), we notice that the potential <fr + evolves in time as A m /o, up to a time-dependent factor F which 
for viable models is practically equal to 1. Such behavior is analogous to that in ACDM when F ~ 1. Of course, in 
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FIG. 5: The parameters Q = (\& — <&) / <3> (left) and 77 = (right) as functions of scale k and redshift z for an /(-R) model 

with w e g = — 1 and = — 10~ 4 . The long-dashed line corresponds to the Compton transition scale Q = 1 (38). One can 
notice the time- and scale-dependent pattern. The parameter uj evolves from u> ~ on scales above Ac (14), to Q ~ 1 on scales 
A < Ac. Analogously, 77 evolves from 77 ~ 1/2 for A 3> Ac to 77 ~ 1 for scales below Ac- 



the case of f(R) the evolution of A m is modified, yet it is an interesting feature of f(R) theories that this relation 
between $ + and A rn is approximately preserved (up to a factor F ~ 1) on sub-horizon scales. A similar result also 
holds for the Modified Source Gravity model [45]. In the Dvali-Gabadadze-Porrati model [46], $ + evolves exactly like 
A m /a [47] without any time-dependent prefactors. 

The change in the relation between the Newtonian potentials $ and has been the focus of some recent litera- 
ture [28, 29]. In [28], a parameter uj was introduced to parametrize the slip between the potentials. In the case of 
f(R) theories, U) is given by the following expression 



* - $ 



JrbSR 
F3> 



(44) 



Its evolution as a function of redshift and scale is shown in Fig. 5. On sub-horizon scales, we can approximately write 
(44) as 



2Q 



2Q + 3 



(45) 



We notice that uj is in general a model, time and scale dependent parameter; on sub-horizon scales, it evolves from 
w~0on scales A ^> Ac, to Q ~ 1 on A <C Ac- 

In [29], an analogous parameter was introduced, 77 = 3 . On sub- horizon scales, this time, scale and model 

dependent parameter assumes the following expression 



n 



2Q 



4Q 



(46) 



Hence it evolves from 77 ~ 1 on scales above the Compton one, to the value 77 ~ 1/2 on scales smaller than the 
Compton scale, as can be seen in Fig. 5. 



3 Note the different sign convention for our <& compared to that in [29], 
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2. In Einstein frame 

It is instructive to derive the above results working in the Einstein frame. In this frame, the modifications are 
described by a scalar field <j> which couples to matter, as described in Appendix B. The perturbation to the scalar 
field, S(f>, is related to the variable \ via the conformal transformation from the Einstein to the Jordan frame, 
x/F — (3nS(f). The behavior of perturbations for a coupled dark energy field has been studied in details in [48], and 
for the case of Chameleon fields in [49]; here we adapt it to the specific case of f(R) theories. The full linear evolution 
equations for the perturbation to the scalar field and for the CDM density contrast are included in Appendix B. In 
this subsection we focus on the sub-horizon regime. 

We notice from (B12), that the effective potential acting on CDM is W — f3K,5(f)/2, where the tilde indicates the 
Einstein frame quantities as in Appendix B. On sub-horizon scales, the equation for the scalar field (Bll) reduces to 
the non-linear Klein-Gordon equation for a massive field, namely 



3/3 a 2 E m S, 



2 k k 2 + a 2 m 2 . 



(47) 



Given this relation between 5(f) and S m , on sub-horizon scales (aH/k <C 1) the dominant contribution to the Poisson 
equation comes from the CDM, and the effective potential acting on a CDM particle is 



(3nS(f> 



3 5^- 
2 k 2 



E S 



1 + 



1 



k 2 



3k 2 + a 2 m 2 



(48) 



Therefore the effect of the scalar field is to introduce a Yukawa type gravitational potential, with a characteristic scale 
A^ = 27r/m0. The mass of the scalar field is uniquely determined as a function of f(R) 



1 FR-AU + R) 
fRR F 2 



(49) 



Therefore the characteristic radius of the Yukawa type potential is A^ ~ (3/rr) 1 / 2 . Let us map the entire equation 
(48) to the Jordan frame. We use the mapping recipe described in [25]; in particular, the relation between the spatial 
curvature perturbations gives ^ = <I' + /3k<5</>/2 and for the energy density of CDM we have S m = 5 m — 2fiK5(j). Finally, 
the scale factors in the two frames are related via a 2 = e^ Kc ^a 2 = Fa 2 . Applying all this to (48), we reproduce eq.(42). 
In particular we see that the physical scale A^, maps to the Compton scale A^ — > (Sfim/F) 1 / 2 w Ac, and that the 
value of the Newton's constant is rescaled as in (43). 



3. Effective dark fluid interpretation 

While working in the Jordan frame, we can consider all the modifications to the Einstein tensor G^ v as components 
of the effective energy-momentum tensor Tff (7). The fluid described by this tensor will not be an independent 
component, but a function of the metric. We can read off the linearly perturbed components of directly from the 
Einstein equations in Appendix A. Given the covariant conservation of the effective energy-momentum tensor (7), we 
need just two quantities to describe the fluid perturbations. We will choose the energy density contrast in the rest 
frame of the fluid, A c g-, and either the anisotropic stress, II ff , or the sound speed, c 2 s , depending on the context. 

The anisotropic stress associated to the effective dark fluid is 



t? tt iu \ 2 k 2 f RR 5R 

E cS U cS (k,a) = - — — - 

9 or b 



2Px 
9 a 2 F 



(50) 



while the effective comoving density perturbation is 

2H 2 



-^^j/rrSR 



(51) 



Equivalently, we can derive the sound speed, defined as 



Ce 2 ff 



5P r 



(c) 



off 



PeffAeff 



(52) 
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FIG. 6: The ratio E e g A c g /(E m A m ) as a function of k and z. At early times and large scales the ratio is very small. After a 
given fc-mode crosses into the Compton radius, the ratio monotonically evolves towards —1/3. The dash line corresponds to 
Q = 1 (38). 



We can now interpret the results of the previous subsections in terms of this effective fluid. On sub-horizon scales 
(k/a ^> H), the effective dark fluid is characterized by the following density perturbation 



Ecfi A c ff 



lfc 2 Q(3F-2) + (F-l) 



3 a 2 



and sound speed 



1 + 3Q 



Q 



(* + *) 



3{l-F)(l + 2Q) + QF 



(53) 



(54) 



The effective shear (50) is proportional to x, therefore its behavior on sub- horizon scales is well described by (39). 
Above the scalaron Compton wavelength, the sound speed is small and positive, c 2 ff ~ + , i.e. the effective fluid 
behaves as an extra sub-dominant clustering component with negligible shear. Eventually, the effective shear II g ~ x 
starts growing and tends to a value comparable to <£>+ on scales well below Ac- The sound speed tends to the negative 
value c 2 ff ~ — 2/[3(3F — 2)] ~ — 2/3 for A -C Ac- A negative value for the sound speed raises concerns about the 
stability of perturbations in the dark component. It would certainly lead to instabilities if there was no anisotropic 
stress associated with the effective fluid. However, in an imperfect fluid with anisotropic stress, one can achieve 
stability even with a negative value of the sound speed. This can be seen easily by combining the evolution equations 
(A1)-(A2) for a generic fluid into a second order differential equation for S. The equation reads 



5" + 3(1 + w)$" 



(1 - 3w) + 



1 



H' 
H 



3(c 2 - w) 1 - 3w 



5' + 3w'$' 




fc 2 / 2 n 



5 + 



i 2 H 2 



(1 + w)y = 



(55) 



The ratio of the anisotropic stress to the energy-density perturbation can be thought of as an anisotropic sound speed, 
c 2 n = II/<5, which, combined with the usual sound speed, can stabilize the perturbations 4 . The stability condition in 



4 It is important to remember that the ratio of the pressure perturbation to the energy density contrast does not necessarily represent 
the physical velocity with which perturbations in the fluid propagate. Perhaps the term isotropic and anisotropic "stiffness" [50] might 
be more appropriate, but we choose to employ the more commonly used term "sound speed". 
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the presence of shear is then generalized to [28] 




(56) 



which is satisfied by the effective dark fluid. 

In Fig. 6 we plot the ratio of the effective dark energy density perturbation to that of the CDM as a function of k 
and z. At early times and large scales, the ratio is very small and A ff is oscillating around zero with a negligibly small 
amplitude. As soon as a given fc-mode crosses into the Compton radius, the ratio begins to evolve monotonically 
towards —1/3. This asymptotic value can also be obtained analytically, by combining eq.(36) with eq.(53) (with 
$+ = ($ + *)/2) to get 



In the limit of Q 3> 1, corresponding to scales well within the Compton radius, the ratio (57) tends indeed to —1/3. 
This value is expected because it corresponds to the change of G e s from G/F to 4G/(3F) (43). 



We have studied the growth of structures in viable f(R) theories. In order to construct f{R) models with specific 
expansion histories we have used the "designer" procedure introduced in [8], and generalized it to include radiation 
and a varying effective equation of state w e s(a). We have shown that viable models, which satisfy all of the conditions 
reviewed in Subsection II B, must closely mimic the ACDM expansion history. In fact, if any departure from w c ff — 
— 1 was observed, it would automatically rule out f(R) as a fundamental theory aimed at explaining the cosmic 
acceleration. One could, in principle, have a viable f{R) cosmology with w c g ^ —1 which would not, however, satisfy 
any of the local tests of gravity. 

The degeneracy with ACDM is broken when one considers the linear perturbations. f(R) models predict a charac- 
teristic scale-dependent growth of large scale structures which may be observationally detectable. To illustrate this 
pattern, throughout the paper we employed the f(R) model with w c s — —1 and = — 1CP 4 . Models with this value 
of would likely fail the solar and galactic constraints which, according to the analysis in [11], require |/^| < 10~ 6 . 
However, the model we used allowed us to demonstrate the qualitative features of the growth pattern better. The 
pattern in the model with = — 1CP 6 would be similar, but with the modifications becoming graphically visible at 
correspondingly smaller scales. Observational signatures of models with |/^| < 1CP 6 would be quite subtle, at a level 
of a few percent difference in the growth factor [11]. Nevertheless, such a level of accuracy should still be within the 
reach of future weak lensing surveys. 

We have solved the exact equations for linear perturbations in the Jordan frame numerically and then confirmed 
our results analytically in the Jordan and Einstein frames using the sub-horizon approximation. We observe a 
characteristic scale-dependent pattern in the growth of structures. The Compton wavelength of the scalaron (14) 
introduces a scale which separates two regimes of sub-horizon gravitational dynamics during which gravity behaves 
differently. On scales A ^ Ac, the scalaron is massive and the "fifth force" is exponentially suppressed, thus deviations 
from GR are negligible. However, on scales inside the Compton radius, the scalaron is light and deviations are 
significant. The relations between $ and *f> , and the relation between them and the matter density contrast, will be 
different below the Compton scale and that affects the growth rate of structures. In particular, on scales below Ac, the 
modifications introduce a slip between these potentials, leading to * ~ 2$. The rate of growth of structures depends 
on the balance between the "fifth force" and the background acceleration. For modes that cross below Ac during 
matter domination the effect of modifications is maximized as the potential $ + grows in the absence of background 
acceleration. When, however, the background starts accelerating, the potentials begin to decay but at a lesser rate 
than in the ACDM model. Therefore, a characteristic signature of f(R) theories would be a non-zero ISW effect 
during the matter era. One should note, however, that the ISW measurements are plagued by large statistical errors 
(because the ISW effect is only a part of the total CMB anisotropy) and cannot provide the percent level accuracy 
needed to test viable f(R) models. Weak lensing studies, on the other hand, are counted on to eventually provide 
highly accurate 3D maps of the gravitational potential $ + . While a comprehensive quantitative analysis of the level 
at which they will be able to constrain is still lacking, our preliminary results, and those of [11, 51, 52] indicate 
that viable f(R) models can still lead to observable differences from ACDM. 

To conclude, we note that the f(R) pattern of structure growth would be similar to that in Chameleon models [49] 
with a universal coupling of the chameleon to all matter fields. These models are somewhat less severely constrained 
than the f(R) models because of the additional freedom in choosing the value of the coupling. 




Q(3F-2) + (F-l) 
1 + 3Q 



(57) 



IV. OBSERVATIONAL OUTLOOK AND CONCLUSIONS 
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Appendix A: Linear Equations in Newtonian Gauge 



The evolution equations for a generic matter field are 



k 



5' + —V - 3(1 + w) — & + 3 — -w)S = 



aH 

V + (1 - 3w)V 



aH 



SP 



5P n 
aH \Sp 5 



bp 



b -(1 + w)* = 

aH 



(Al) 
(A2) 



The linear Einstein equations are 
• time-time 



1,2 op f u2 \ 

6* + 6$' + = -—±Si - f R ( 6* + 6$' + 2^2$J + 3f RR bR' - f' R (6* + 3$') + 



3f RR 1 



W_ 
~~H 



SR ~ ^2 H ^fRR + 3 fRR 



5R 



momentum 



space- diagonal 
H' 



k q, + = - f R (fc* + fc$') + ±k(f RR 5R)' - \kf RR 5R - ^kf' R * 



(A3) 



(A4) 



3 + 2' 



H 



1 k 2 



3a 2 H 2 



* + + ( 1 + 4" ) + 2*' + + 



H 



1 k 2 



3Ei SR 



3a 2 H 2 2H 2 bp 



5 + 



3 + 2 



~H 



lk 2 .. ( H' 

# - -^TTTTT^ + $ + 1 



3a 2 H 2 



H 



1 k 2 

$' + 2$' + ii'' + -^_ T $ 



3a 2 H 2 



- 3 



H 



+ 



f RR {5R)" 



+ 1 



+ 



+ 



3a 2 ff 2 



H 



-^(^' + 2* + 2$') 



• space off- diagonal 



$ - * = Ip-BiHi - /r ($-*) + /*r£R • 



The linear perturbation to the Ricci scalar is 

k 2 



5R 



-2H Z 



k 2 

2 $ 
a 2 # 2 a 2 # 2 



# + 3$" + W + 3$' + 6 2 + — * + 3 1 + — W>' 



H 



H' 



H 



(A5) 



(A6) 



(A7) 



Appendix B: Mapping to the Einstein Frame 

Using the approaches of Chiba [20] and of Magnano & Sokolowski [19], following [30], we recast the gravitational 
action (1) into a dynamically equivalent form by introducing an intermediate scalar field The equivalent action is 
[19] 

S =J^J d 4 x^—g [(<& + /(<&)) + (1 + U)(R - *)] + J d^x^—gC^g^v] , (Bl) 
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where /$ = df/d&. One can verify that, if d 2 f/d& 2 7^ 0, the field equation for $ is R = <f>, which reduces (Bl) to 
the original action. Next consider the conformal transformation 



such that the function oj(x a ) satisfies 



e- 2 "(l + f R ) = l 



(B2) 



(B3) 



With this choice of uj the action (Bl) transforms into an action with the usual Hilbert-Einstein form for gravity. If 
we now define the scalar field <fi = 2u>/0k, where /3 = y/2/3, the resulting action becomes 



S = ^ J d A x^gR + J d A x^g 



+ J d A x^e~ 2 ^C m [^e^g^} 



where the potential V(<j>) is determined entirely by the original form (1) of the action and is given by 

v{<t,) = 2^oTW 

The Einstein-frame line element can be written in familiar FRW form as 

dS 2 = d 2 (-dr 2 + dx 2 ) , 



(B4) 



(B5) 



(B6) 



where the Jordan and Einstein metrics are related through the conformal transformation a 2 = e^^a 2 . It is also 
convenient to define an Einstein-frame matter energy-momentum tensor by 



f^ = {~p + P)U^U v + P~g^ , 

where % = e^/ 2 U^p = e^^p and P = e~ 2 ^P. 

The equations of motion obtained by varying the action with respect to the metric g^ v are 



(B7) 



(B8) 



and arc more familiar than those in the Jordan frame, although there are some crucial distinctions. Most notably, 
in this frame test matter particles do not freely fall along geodesies of the metric g^ v , since the scalar field is also 
coupled to matter. 

The remaining equations of motion, for the scalar field and for the perfect fluid matter, are given respectively by 



<j) + 2H4> + 5% = -nfra 2 {p - 3P) , 

p + m{p + P) = ~K^{p-w) . 

where = dVj d(f>. 

At linear level in the perturbations, the evolution equations for scalar field and CDM perturbations are 



61 + 2 + ^ 



+ 



~a 2 H 2 



H 2 



H 



3/3 ft j 

- X —On 

2kH 2 

+ T 



= 



The mapping prescription at the linear level is the following 
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